Electron spin resonance detected by a superconducting qubit 
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A new method for detecting the magnetic resonance of electronic spins at low temperature is 
demonstrated. It consists in measuring the signal emitted by the spins with a superconducting 
qubit that acts as a single-microwave-photon detector, resulting in an enhanced sensitivity. We 
implement this new type of electron-spin resonance spectroscopy using a hybrid quantum circuit in 
which a transmon qubit is coupled to a spin ensemble consisting of NV centers in diamond. With 
this setup we measure the NV center absorption spectrum at 30 mK at an excitation level of ~ 15 \±b 
out of an ensemble of 10 11 spins. 



I. INTRODUCTION 

Electron spin resonance (ESR) spectroscopy at low 
temperatures is often complicated by the long spin-lattice 
energy relaxation time which can reach minutes at sub- 
kelvin temperatures. For continuous- wave (CW) ESR 
spectroscopy, this implies working at low powers to avoid 
saturating the spins; in pulsed ESR this imposes low rep- 
etition rates. In both cases a higher sensitivity for detect- 
ing the signal absorbed or emitted by the spins would 
be desirable. Recently, tools borrowed from supercon- 
ducting quantum electronics have been applied to high- 
sensitivity ESR spectroscopy. High-Q superconducting 
coplanar waveguide resonators have been used for CW- 
ESR at millikelvin temperatures with a cryogenic low- 
noise HEMT amplifier followed by homodyne detection 
[1-4J, and for pulsed ESR at kelvin temperatures with a 
commercial ESR spectrometer |5J , yielding promising re- 
sults in terms of sensitivity. Here we go one step further 
and use an on-chip single microwave photon detector [5] 
based on a superconducting qubit [7] to realize a new type 
of high-sensitivity low-temperature ESR spectrometer. 

The principle of our experiment is compared to more 
conventional ESR techniques in Fig. [T] In conventional 
ESR spectroscopy, a microwave pulse is applied to an 
ensemble of spins close to their resonance frequency u s 
through a low-Q cavity of frequency uj c with which the 
spins are tuned in resonance (uj s = uo c ) by a magnetic 
field. After being excited by the pulse, the spins re- 
radiate coherently part of the absorbed energy through 
the cavity into the detection waveguide, giving rise to a 
free induction decay (FID) signal measured by homodyne 
detection, which yields after Fourier transform the spin 
absorption spectrum [8j . What limits the sensitivity of a 
typical commercial CW-spectrometer operating at 300 K 
to ~ 10 10 spins/V / Hz for a linewidth of 0.1 mT and an 
integration time of 1 s is the overall noise temperature 
of the detection chain. In this work, we replace the de- 
tection chain by a superconducting qubit and its readout 




Figure 1. Comparison between conventional pulsed electron 
spin resonance (ESR) and qubit-detected ESR. a, Conven- 
tional ESR: A spin ensemble is placed in a cavity with fre- 
quency ujc and driven with a microwave pulse resonant with 
the cavity. When the ESR frequency uj s matches cj c , the 
spins absorb the microwave pulse, and emit immediately af- 
ter a free-induction decay (FID) signal into the waveguide 
connected to the cavity, b, Qubit-detected ESR: The cavity 
is now frequency-tunable and embeds both the spin ensemble 
and a superconducting qubit with frequency uj ge . In a first 
step (1), the spins are probed by a spectroscopy pulse which 
excites them if its frequency matches uj s . In a second step (2), 
the cavity frequency is tuned to uj c — oo s , receives the FID sig- 
nal from the spins, and is afterwards tuned to transfer this 
signal to the superconducting qubit at uj ge . Finally (3) the 
qubit excited state probability P e (cj s ) is measured, mapping 
the spins absorption spectrum. Very low excitation powers 
can be used given the high sensitivity of the method. 



circuitry. This results in an increased sensitivity since a 
superconducting qubit is a nearly ideal single microwave 
photon detector |6J at its resonance frequency uo ge . In 
order to transfer part of the excitation of the spins to 
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Figure 2. a, Sketch of the implementation of qubit-based 
ESR. The spin ensemble S consists of NV centers in a di- 
amond crystal. They are coupled to the frequency-tunable 
coplanar waveguide resonator C used as the ESR cavity. C 
also embeds the ESR-detector qubit Q, a superconducting 
qubit of the transmon type whose state can be readout with 
another resonator R. Microwave pulses for spin spectroscopy 
as well as for qubit readout are sent via an input port cou- 
pled both to C and to R. b and c, Sketch and energy levels 
of NV centers in diamond. In our setup a DC magnetic field 
Bjvy is applied along the [1,1,1] direction, resulting in dif- 
ferent Zeeman splittings for centers having the N — V axis 
parallel T$nv (ensemble /, in red) and those having their axis 
along the three other (1,1,1) axes (ensemble III, in blue). 
The ESR frequencies uj±iju are further split in three reso- 
nance lines due to the hyperfine interaction with the spin-1 
14 N nuclear spin |11| . 



the superconducting qubit, the ESR resonator is made 
frequency-tunable and with a high quality factor. Spec- 
troscopy is performed by first exciting the spins with a 
weak microwave pulse, collecting the radiated FID signal 
with the resonator tuned at uo Sl then transferring this 
signal to the qubit at uo gei and measuring its final state. 
Repeating this experimental sequence yields the proba- 
bility P e to find the qubit in its excited state, which re- 
produces the spin absorption spectrum. The sensitivity 
of such an ESR spectrometer is set by the efficiency at 
which signal photons can be transferred from the spins 
to the resonator, then to the qubit, and by the fidelity 
with which the qubit state can be measured. 



II. DEVICE AND EXPERIMENTAL SETUP 

We implement this method using a recently reported 
hybrid quantum circuit [9J that includes an ensemble of 
electronic spins magnetically coupled to a superconduct- 
ing resonator, itself electrically coupled to a supercon- 
ducting qubit, as sketched in Fig. [2ji. The spins are 
negatively-charged nitrogen-vacancy (NV) color centers 



in diamond, whose structure and energy levels are sum- 
marized in Fig(2]3 and c. The ground state of NV centers 
has a spin one with splitting uj±/2tt ~ 2.88 GHz between 
states ms = and ms — ±1 at zero magnetic field [10J. 
Each of the two ms = to ms = ±1 transitions is further 
splitted into three peaks separated by 2.2 MHz due to the 
hyperfine (HF) coupling to the 14 N nuclear spin jTT] . In 
the experiment a static magnetic field Bnv = 1.1 mT [12] 
is applied along the [1,1,1] crystallographic axis to lift 
the degeneracy between states ms = ±1. Centers having 
their N — V axis along [1, 1, 1] (called in the following 
ensemble I) undergo a different Zeeman shift from those 
along the three other (1, 1, 1) axes (ensemble ///), result- 
ing in two different ESR frequencies cj + i/27t = 2.91 GHz 
and co> + ni/27r = 2.89 GHz for the ms = to ms = +1 
transition on which we will exclusively focus in the fol- 
lowing. 

The diamond crystal used is of the High-Pressure High- 
Temperature (HPHT) lb type and has a NV center con- 
centration of ~ 3 ppm and a residual nitrogen concen- 
tration of ~ 30 ppm. In our setup, it is glued on top of 
the ESR cavity C, a coplanar waveguide superconduct- 
ing resonator [13J of quality factor Q ~ 10 4 made of a 
niobium thin-film sputtered on a silicon substrate. The 
spin ensemble S detected in the experiment consists of 
the ~ 10 11 NV centers that lie within the mode volume 
of C, thus within a few microns of the diamond surface. 
The cavity frequency uj c (&) can be tuned on a nanosec- 
ond timescale by a flux $ applied through the loop of a 
superconducting quantum interference device (SQUID) 
inserted into C [U fl4HT6] . The superconducting qubit 
is a Cooper-pair box of the transmon type [T7l [T8] with 
resonance frequency uj ge between its ground state \g) and 
excited state \e). It is coupled to an additional resonator 
R which is nonlinear and used to read-out the qubit state. 
As explained in details in Ref. [19], this readout is per- 
formed in a single shot by measuring the phase of a mi- 
crowave pulse reflected on R. This phase takes two dif- 
ferent values depending on the qubit state, and repeating 
~ 10 4 times the same experimental sequence yields the 
qubit excited state probability P e . 



III. SINGLE PHOTON STORAGE AND 
RETRIEVAL 

As reported in an earlier work |9J (see also [21J), it is 
possible with this circuit to coherently exchange a sin- 
gle quantum of excitation between the spin ensemble S 
and the qubit Q via the tunable cavity C. To demon- 
strate that, the qubit is prepared in |e); its excitation is 
transferred to the cavity by sweeping adiabatically ou c (&) 
through uo ge , which is then tuned suddently in resonance 
with the spins at ujk (where K = + J, +11 1) for some 
interaction time r. The excitation left in the cavity is fi- 
nally transferred back into the qubit, which is then read- 
out. As shown in Fig. [3j the resulting qubit excited 
state probability P e {j) is found to oscillate, revealing the 
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Figure 3. Single photon transfer between qubit and spins 
[9]. (inset): Pulse sequence used to excite the qubit in state 
\e) and transfer its excitation to S via C back and forth, 
(panel): Swap oscillations for the two spin frequencies cj+i 
and Co>+m . Swap times r s ,i and r s ,in are indicated by arrows. 
Note that the two curves are shifted vertically for clarity with 
corresponding scales on left (for ensemble 7) and right axes 
(ensemble 777). 



conversion of a single microwave photon into an elemen- 
tary collective excitation of the spin ensemble. For well- 
defined interaction times t s ^k (see Fig. [3|, the excitation 
in the qubit is swapped into the spin ensemble |9J; at 
a later time it is recovered in the qubit with a fidelity 
- 0.1. 

A quantum-mechanical description of this experiment 
[22, 23J is useful in the discussion of the ESR results 
presented below. Each of the Nk effective spins of the 
ensemble at ujk is modelled as an effective harmonic os- 
cillator with frequency Uj K and annihilation (resp. cre- 
ation) operator bj K (resp. bj K ), an approximation valid 
in the low-excitation limit as is the case throughout 
this article. The spin ensemble and cavity are then de- 
scribed by Hamiltonians ^2hujj K b^ K bj K and huj c (<fr)a^a, 
a (resp. a)) being the cavity annihilation (resp. cre- 
ation) operator. The coupling between the resonator 
and the spin ensemble is described by a Hamiltonian 
Hk = — 9j K (^] K a + h.c), with gj K the coupling 
constant of spin jx to the resonator. This Hamilto- 
nian can be rewritten as Hk = —i^gx (^ K a + h.c?j with 

1/2 

qk = (J2d] K ) the spin ensemble-resonator collective 
coupling constant and bx = (1/gx) X^jK^K ^ ne annihi- 
lation operator of the collective spin excitation coupled to 
the cavity. This superradiant mode has a spatial profile 
given by the spatial dependence of the coefficients 
which reproduces the profile of the magnetic field inside 
the cavity mode. Note also that bx involves all the spins 
belonging to group K, even if they have different fre- 
quencies due to slightly different magnetic environment 
in the crystal, also including the three possible states of 
the 14 N nuclear spin causing the hyperfine structure; as 



a result this mode is coupled to Nk — 1 dark modes that 
act as a bath [16j |22 -24J. Using these notations, one 
describes the oscillations shown in Fig. [3] as occurring 
between states |1 c ,0k) and |0 C , 1^), where |l c ) = a) |0 C ) 
is the usual Fock state with 1 photon in the cavity, and 
I Ik) = | Ox) is the first excited state of the superradi- 
ant mode; damping of these oscillations is due to inhomo- 
geneous broadening, and can be interpreted as damping 
of state \1k) into the bath of dark states [IB] . 



IV. ESR PROTOCOL AND DISCUSSION 

The ESR protocol is shown in Fig. [4] It consists in an 
experimental sequence similar to the one used for single 
photon storage (see fig. [3|, but with the spin ensem- 
ble initially excited at several photons level instead of 
the qubit prepared in its excited state. More precisely, 
a low-power microwave pulse of duration At = 2 jis and 
varying frequency oo p is applied to the spins while uj c (&) 
is far detuned; the resulting excitation is transferred first 
into the cavity by tuning uj c (&) suddenly in resonance 
with ojk for the swap time t s ^k, then into the qubit by 
an adiabatic swap interaction; the qubit state is finally 
measured. Provided the average number of microwave 
photons emitted by the spins into the ESR cavity stays 
much lower than 1 to avoid saturating the qubit, the re- 
sulting excited state probability P e (uo p ) is expected to 
reproduce the spin ensemble absorption spectrum. Ex- 
perimental results of Fig. [4Jd indeed display the charac- 
teristic HF structure of NV centers consisting in three 
peaks separated by 2.2 MHz for both spin ensembles +7 
and +777. This validates the concept of electron spin res- 
onance detected by a superconducting qubit. Note that 
P e (u p ) <C 1 showing that qubit saturation is avoided as 
wanted. 

We now discuss the sensitivity of this qubit-based ESR 
spectrometer. The qubit state is detected in a single-shot 
with a fidelity of ~ 0.7 at the end of an experimental 
sequence that lasts typically 50 /is, yielding a 1% preci- 
sion on the probability P e in one second. To translate 
this sensitivity in a magnetic moment unit, one needs to 
know with what efficiency the excitation of the spin en- 
semble is actually transferred to the qubit. The transfer 
of one microwave photon from the cavity to the qubit is 
performed with an efficiency of order unity (in our ex- 
periment it is around 0.7 limited by losses in the cavity 
and qubit), so the limiting factor is the efficiency of the 
transfer of the spin ensemble excitation to the cavity dur- 
ing their resonant interaction. At first sight one might 
think that since the spin ensemble and cavity are in the 
strong coupling limit, one excitation of the spin ensemble 
should also be converted into a microwave photon with 
an efficiency of order 1, similar to what happened in the 
coherent oscillations shown in Fig. [3] This reasoning 
is not correct here because the collective spin mode b UJp 
excited by the spectroscopy pulse does not necessarily 
match perfectly the superradiant mode bK- Indeed, al- 
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Figure 4. a, Experimental pulse sequence used for qub it- 
detected ESR: (1) The spins are first weakly excited by a 2 /is 
microwave pulse with a frequency uj p ] (2) the resulting spin 
excitation is transferred to the cavity C by a fast flux pulse 
which brings u c in resonance with ujk {K — +1, +11 1) for 
a swap time t s ,k, and then to the qub it Q by an adiabatic 
swap (aSWAP). (3) The qubit excited state probability P e 
is finally measured, b, Measured (open circles) and calcu- 
lated (solid line) P e (uo p ) for spin ensemble +11 1 (left) and 
+/ (right). The spin density p(uj) used in the calculation is 
shown as a dashed line. 



though the spatial matching of the two modes is excellent 
since the specroscopy pulse is applied through the cavity, 
this is not the case for spectral matching: only spins hav- 
ing a resonance frequency within the spectroscopy pulse 
bandwidth S/2tt = 150 kHz around uo p contribute to 6 W , 
whereas all spins within the hyperfine line (total width 
A/27T ~ 5 MHz) contribute to the superradiant mode bx- 
As a result, one expects an overlap of order y/5/A be- 
tween the b 0Jp and bx modes, implying that there should 
be ~ A/S = 20 times less excitations in the bx mode 
(and thus also in the cavity mode after the swap interac- 
tion) than in the b 0Jp mode. 

This argument can be made rigorous and 
quantitative. The collective spin mode ex- 

cited by the spectroscopy pulse is defined as 

K = [T,9jK a u; p (u jK )b jK ] /y / E#L K P (^)| 2 > 

with a Up (uj) the pulse Fourier transform. The 
quantity of interest is then the correlation function 

f3(cjp,T 3 ) — (a(T 3 )blj p (0)^ giving the probability ampli- 
tude for an excitation created in b^ by the spectroscopy 



pulse to be transferred into a photon inside the cavity 
mode after an interaction time r s . This function can be 
computed numerically given a certain spin distribution 
p(uj) using the formulas derived in the Appendix. 
In our experiment, the linewidth of each hyperfine 
peak w +i /2tt = 1.6 MHz and w +iii /2tv = 2.4 MHz, 
and coupling constants g+i/27r = 2.9 MHz and 
g+iu/27r = 3.8 MHz have been determined from other 
measurements [9j, so that a direct comparison with 



theory without any adjustable parameter is possible as 
shown in Fig. [4j The agreement is quantitative (note 
that we have also included in the distributions p(uj) 
additional ESR frequencies caused by the hyperfine 
interaction of the NV center with neighboring 13 C nuclei 
with the 1.1% natural abundance as expected). From 
this calculation, we deduce that the average excitation 
of the spin ensemble at resonance in the data shown 
in Fig. [4] is ~ 15, in agreement with the qualitative 
argument presented above. In the present state of the 
experiment, the qubit-based ESR spectrometer therefore 
measures the spectrum of an ensemble of 10 11 NV 
centers at an excitation level of order 15 in a one 
minute integration time. 

Thanks to this very low excitation level, the exper- 
imental sequence can be safely repeated at 20 kHz, de- 
spite the NV centers energy relaxation time reaching min- 
utes at 30 niK [3l [20] . More precisely, two factors con- 
tribute to make this experiment possible: 1) at the end 
of each experimental sequence the excitation of the b 0Jp 
mode quickly decays into the bath of dark modes, allow- 
ing the next experimental sequence to start with b UJp in its 
ground state and thus keeping the average number of ex- 
citations transfered to the qubit well below 1 as needed 
to avoid saturation, 2) the low excitation rate ensures 
on the other hand that the ensemble of 10 11 spins stays 
far from saturation even after repeating the sequence for 
hours. 

We finally note that our calculation reproduces a puz- 
zling feature of the data that was not discussed yet: 
the middle peak of the P e {oo p ) curve has a lower am- 
plitude than the two other peaks, both for the +1 and 
the +III curves as seen in Fig. [4j although the spin den- 
sity p(uS) used in the calculation is a simple sum of three 
Lorentzians with the same amplitude. Our ESR protocol 
thus appears to slightly distort the absorption spectrum. 
This phenomenon originates from the uo p dependence of 
the energy transfer efficiency from the spins into the cav- 
ity, caused by the fact that gx ~ A in our sample. It 
could probably be corrected in future experiments either 
by increasing gx or by transferring the spin excitation to 
the cavity with an adiabatic passage. 

Besides detecting a large ensemble of N = 10 11 elec- 
tronic spins at near single-excitation level, it is interest- 
ing to discuss what is the minimal number of spins N m i n 
that could be detected with a similar experimental pro- 
tocol in order to compare it to the sensitivity of existing 
conventional spectrometers. For that we will change per- 
spective in the following discussion, and assume that the 
spins being measured can actually be excited at satura- 
tion. We suppose that the N spins, of inhomogeneous 
linewidth A, have been excited by a hard ir/2 pulse. In 
the weak coupling limit g\/N <C k, <C A (k, = uo c /Q be- 
ing the cavity damping rate), the spins emit in the cavity 
ft = g 2 N 2 /(4A 2 ) photons (see Appendix B). Taking a 
conservative estimate for the minimal average excitation 
that can be detected by a superconducting qubit within 
one second to be 0.05, a spin-cavity coupling constant 
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g/2n = 10 Hz, one obtains N min = 10 5 spins/A/Hz for 
a 0.1 mT linewidth corresponding to A/2tt = 2.8 MHz. 
This figure is five orders of magnitude better than a com- 
mercial spectrometer at 300 K, and two orders of magni- 
tude better than the record sensitivity of 10 6 spins/ VHz 
for a 0.01 mT linewidth that was recently reported with 
a surface loop-gap resonator operated at 10 K [26] and a 
coplanar waveguide resonator at 4 K |5J. Note however 
that in order to operate such a qubit-based spectrometer 
in practice, one would need 1) to use a repetition rate 
around 10 kHz which requires a reasonably short spin- 
lattive relaxation time or some way to repump rapidly 
the spins into their ground state (optically [25] or elec- 
trically for instance) and 2) to find a qubit design that 
withstands large magnetic fields usually needed for ESR. 
Interesting alternative possibilities could be to physically 
separate the qubit-detector from the spins, which would 
allow more easily to apply large magnetic fields to the 
spins without perturbing the qubit, and to continuously 
monitor the qubit state with a parametric amplifier as 
demonstrated in recent experiments [27] instead of puls- 
ing the qubit state detection as done here. 

In conclusion we have discussed a new type of ESR 
spectrometer in which the signal coming from the spins 
is detected by a superconducting qubit acting as a single- 
microwave-photon detector. We have implemented this 
idea on an ensemble of ~ 10 11 NV centers coupled to 
a transmon qubit, measuring their absorption spectrum 
at an excitation level of ~ 15 with a well-resolved 
hyperfine structure. Estimates indicate that this spec- 
trometer would be able to detect 10 5 spins/ a/Hz with a 
0.1 mT linewidth, a gain of two orders of magnitude in 
sensitivity compared to the best reported values for con- 
ventional spectrometers. Our work thus demonstrates 
the potential of superconducting circuits for electron spin 
resonance spectroscopy. 



APPENDIX A 

We now explain in more detail how the theory curves 
in figure 3 are calculated. As explained in the main 
text, the energy transfer efficiency from the resonator 
to the qubit is of order unity and can be well mod- 
eled by an ideal adiabatic passage. In this context 
the quantity of interest is the resonator population af- 
ter the interaction with the spins. The calculations 
are performed in the Holstein-Primakoff approximation, 
in which the spins and the resonator are described by 
harmonic oscillators. The total system Hamiltonian is 
H/h = u c (<$>)a^a + E^j^j + T,i>9i0>]a - b^), g } be- 
ing the coupling constant of spin j with the resonator. 
We need to calculate the probability that the excitation 
created at t = in the spins to be transferred to the 
cavity after a time t, this probability is the square mod- 
ulus of (0| a(t)bl Jp |0). The spins excitation is created by 
a microwave pulse of central frequency uj p with a pulse 
envelope in frequency described by a Up (uo) = a(uj — oj p ). 



a typical envelope is a Lorentzian function with FWHM 
6. We can define an operator ftj^ that describes the ex- 
citation induced by this pulse as 
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this comes simply from the standard atom-field inter- 
action for a classical light source such as the one used 
in the experiment. As shown in [22J the quantity 
(0| |0) can be calculated by considering an effec- 

tive non-Hermitian Hamiltonian 



Heff/fr- 



( 



V 



w igi igi 
-igi tii 

-ig 2 u> 2 



\ 



■J 



(2) 



with complex angular frequencies ujq = uj c (<&) — Ik/2 and 
ujk — ojk — ijo/2 ; here, 70 is the spontaneous emission 
rate of each spin (that we take here to be zero since NV 
centers at low temperature have negligible energy relax- 
ation). Indeed, introducing the vector X(t) of coordi- 
nates \(a(t) (0) ),..., (bj(t)aJ (0)) , ...] it can be shown 
that dX/dt = —(i/H)H e ffX. The formal solution to this 
equation is then 



X(t) = C- 1 [(s + iH eff /h)- 1 X(0)], 



(3) 



which gives <0| a(t)bt p |0) = x G * • X(t) = C' 1 [t Up (s)] 
with x G = (1, 0, 0, ...) and C[f(s)] = f e~ st f(t)dt (s being 
a complex number). The initial condition X(0) is the one 
produced by b^ given in EqJTJ thus 
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T,k a ^k)g 2 k 



i 7 o + ^/EjKM)\ 2 9j ' 

where t\{—iuj) = ij [cj — ojq + w/2 — W(lj)] with 
W(uj) = Ej 9j/ i 00 ~ 00 j + 270/2]. Note that we evalu- 
ated tcj (s) for s = — iu, this is sufficient to perform the 
Laplace transform inversion as there are no singularities 
in the imaginary axis of t 0Jp . We define the spin density 

p(uj) encompassing the coupling strength, which is pos- 

2 

sibly different for each spin, as p(uj) = E 7 - S(u — ujj). 
Using this definition in the equation above we have 



g K h(-iuj) (a*p)(u p ) 



ijo + - Up) yj(\a\ 2 *p)(u p ) 



(5) 
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The spectral width of the microwave pulse is, in our 
case, much smaller than any scale that characterizes our 
distribution p(ui). This allows the rewriting of the con- 
volution above as 



(a*p)(u p ) 

\/(M 2 *p)(up) 



(6) 



where the constant A 



f ot(uj)duj 



is purely character- 



ized by the pulse envelope with no dependence on u p , 
yielding for example A = V^\AV^ f° r a Lorentzian en- 
velope. This means that if we consider that the spins 
are distributed at a typical range A the equation above 
gives a rigorous justification of the rule of thumb that 
says that the efficiency of the spins-resonator transfer is 
given by the overlap ^/ 5 / A. 

Finally to generate the theoretical curve in Fig. 3, we 
perform a numerical inversion of the Laplace transform 



for each uj p and take 
t = t S j. 



<0|a(t)6t | )| 2 at t 



ts,iii or 



APPENDIX B 

We now explicit the calculation of the sensitivity of 
our qubit-based ESR spectrometer in the weak coupling 



limit qk <C k <C A . The Hamiltonian coupling the 
spins to the cavity field is H = hg(S-a) + h.c), where 
S- — (J i,- Dem g the lowering operator of spin 

i. In the absence of driving field, the equation for the 
intra-cavity mean field is then easily obtained as 



d < a> 
dt 



< a > -ig < S- > . 



Right after a tt/2 pulse on the spins, | < S- > | = N/2. 
Neglecting the back-action of the cavity field on the spins 
(which is justified in the weak coupling limit), we simply 
get that < S- >= (N/2)e~ At . From that one shows that 
<a> (t) = -igN(e- Kt l 2 - e~ At )/[K - 2A] which in the 
limit k <C A yields a maximum photon number in the 
cavity of n = g 2 N 2 /(4A 2 ). 

Acknowledgements We acknowledge useful discus- 
sions with B. Julsgaard, K. Moelmer, J. Morton, A. 
Briggs, and within the Quantronics group, and techni- 
cal support from P. Senat, P.-F. Orfila, T. David, J.-C. 
Tack, P. Pari, P. Forget, M. de Combarieu. We acknowl- 
edge support from European project Solid, ANR projects 
Masquelspec and QINVC, C'Nano, Capes, and Fondation 
Nanosciences de Grenoble. 



[1] Y. Kubo et al., Phys. Rev. Lett. 105, 140502 (2010). 
[2] D.I. Schuster, et al., Phys. Rev. Lett. 105, 140501 
(2010). 

[3] R. Amsuss, et al., Phys. Rev. Lett. 107, 060502 (2011). 

[4] P. Bushev et al., Phys. Rev. B 84, 060501 (2011). 

[5] H. Malissa, D. I. Schuster, A.M. Tyryshkin, A. A. Houck, 

and SA. Lyon, arxiv:1202.6305. 
[6] Y.-F. Chen et al., Phys. Rev. Lett. 107, 217401 (2011). 
[7] G. Wendin and V.S. Shumeiko, Low Temperature Physics 

9, 724 (2007). 

[8] Note that in a real ESR spectrometer, the FID signal 
would in fact be difficult to detect using this protocol 
because of ringing of the cavity and of the detectors 
dead-time. The spin absorption spectrum would therefore 
rather be measured by CW-ESR rather than by this FID- 
detected protocol. We nevertheless discuss it for pedagog- 
ical reasons, because it is closest to our qubit-detected 
ESR experiment in which these technical issues are irrel- 
evant. 

[9] Y. Kubo et al., Phys. Rev. Lett. 107, 220501 (2011). 
[10] F. Jelezko, T. Gaebel, I. Popa, A. Gruber, and J. 

Wrachtrup, Phys. Rev. Lett. 92, 076401 (2004). 
[11] S. Felton, et al, Phys. Rev. B 79, 075203 (2009). 
[12] The experiment was performed in the same magnetic field 



as in [9 where the reported value was mistakenly men- 
tioned to be Bnv = 1.4 mT. 
[13] L. Frunzio, A. Wallraff, D. I. Schuster, J. Majer, and R. J. 
Schoelkopf, IEEE Trans, on Applied Superconductivity 
15, 860 (2005). 

[14] A. Palacios-Laloy et al, J. Low Temp. Phys. 151, 1034 
(2008). 

[15] M. Sandberg, et al, Appl. Phys. Lett. 92, 203501 (2008). 

[16] Y. Kubo, et al, Phys. Rev. A 85, 012333 (2012). 

[17] J. Koch et al, Phys. Rev. A 76, 042319 (2007). 

[18] J. A. Schreier et al, Phys. Rev. B 77, 180502 (2008). 

[19] F. Mallet, et al, Nature Phys. 5, 791 (2009). 

[20] A. Jarmola, V. M. Acosta, K. Jensen, S. Chemerisov, D. 

Budker, arxiv:1112.5936 
[21] X. Zhu et al, Nature 478, 221 (2011). 
[22] I. Diniz et al, Phys. Rev. A 84, 063810 (2011). 
[23] Z. Kurucz, J.H. Wesenberg and K. Molmer, Phys. Rev. 

A 83, 053852 (2011). 
[24] K. Sandner et al, arxiv: 11 12.4767 (2011). 
[25] A. Yang et al, PRL 102, 257401 (2009). 
[26] Y. Twig, E. Dikarov, W.D. Hutchison, and A. Blank, 

Rev. of Sci. Instr. 82, 076105 (2011). 
[27] R. Vijay et al, Phys. Rev. Lett. 106, 110502 (2011). 



